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^ ; Recently Shukurov et al [Phys Rev E 72, 025302 (2008)], made use of non- 

^ ' orthogonal curvilinear coordinate system on a dynamo Moebius strip flow, 

^ to investigate the effect of stretching by a turbulent liquid sodium flow. In 

; plasma physics, Chui and Mofl"att [Proc Roy Soc A 451,609,(1995)] (CM), 

■ considered non-orthogonal coordinates to investigate knotted magnetic flux 

tube Riemann metric. Here it is shown that, in the unstretching knotted 
tubes, dynamo action cannot be supported. Turbulence there, is generated 
by suddenly braking of torus rotation. Here, use of CM metric, shows that 
stretching of magnetic knots, by ideal plasmas, may support dynamo action. 
Investigation on the stretching in plasma dynamos, showed that in difl"usive 
media [Phys Plasma 15,122106,(2008)], unstretching unknotted tubes do not 
support fast dynamo action. Non-orthogonal coordinates in flux tubes of non- 
constant circular section, of positive growth rate, leads to tube shrinking to a 
constant value. As tube shrinks, curvature grows enhancing dynamo action. 
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I Introduction 



Earlier R Ricca [1 ] has developed an interesting method of investigating the vortex fila- 
ments by applying a very popular method from Einstein gravity [2], called Ricci rotation 
coefficients to obtain covariant or intrinsic equations equivalent to Da Rios equation and 
applied to solitons in n-dimensional space. The manifold considered is a not necessarily, 
torsion-free Riemannian space endowed with a Levi-Civita connection. In this paper, 
attention is drawn to a simple analytical Riemannian geometrical model, of a metric of 
a knotted magnetic flux tube, previously investigated by Chui and Moffatt [3]. In their 
work they considered applications to magnetic helicity and energy, but did not take into 
account implications of their metric for dynamo action. In the present paper one pretends 
to fill this gap by considered the Ricca's 3D RRCs tool to obtain an invariant classifi- 
cation of knotted magnetic flux tube in its unstretching and non-dynamos aspects, that 
allows us to associate for example, the vanishing of some of the components of RRCs and 
the non-existence of dynamo action. In principle the reasoning here applies to kinematic 
dynamos. Other covariant model important in magnetic reconnection of astrophysical 
plasmas has been recently developed by Titov et al [4]. The ideas of simple dynamo 
models in Riemannian space were inspired by the cat Arnold's dynamo first toy model 
[5] uniformly stretched in Riemannian space. More recently one has investigated the ap- 
plication of Vishik [6] anti-fast dynamo theorem to the dynamo plasma, showing that 
only a slow plasma dynamo is supported when tubes are unstretched. Though most 
of stellarators and tokamaks are metalic unstretched tubes their dynamo magnetic flow 
inside maybe stretch, so in principle fast dynamos may happen in tokamaks even in Eu- 
clidean spaces. Recently Nornberg at al [7] have shown that the stretch-twist and fold 
fast dynamo method [8] could be obtained from Madison torus experiment. The paper is 
organized as follows: In section II the model is presented and the unstretching condition 
is imposed on the magnetic field flow. Section III presents the computations of stretching 
Ricci rotation coefficients and the magnetic energy of helical knotted Riemannian tube 
manifolds. In section IV it is shown, that in the case of orthogonal and non-orthogonal 
coordinates, dynamo action can be present. Conclusions are presented in section V. 
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II Riemann-Chui- Moffat t knotted magnetic flux tube 
metric 



This section presents the Riemann metric proposed by Chui and MofFatt of a curvilinear 
coordinates {s, r, (j)). Note that magnetic lines along the tube are computed in the Frenet 
frame (t, n, b) which obeys the following evolution equations 

| = «M. (II.1) 

^ = -K(s)t + rb (II.2) 
as 

f = (II.3) 

where, t is the vector tangent to magnetic lines while in the plane orthogonal to them, 
where the normal n and the binormal vector b are defined. Here n and r are Frenet 
curvature and torsion scalars.An arbitrary vector x in the tube T can be expressed as 

x{s) = X + r(cos en + sin Ob) {II A) 

Where (r, 9) are polar coordinates that are defined in the plane through X(s) defined by 
the vectors (n, b). The boundaries of a ribbon surface C — 9 — constant along the tube, 
and its curve of intersection Ve with the boundary tube &T . Here M is the Gauss linking 
topological number of (C, Tg) which my be positive, negative or even zero. Now, let us 
define a twist angle 

(11.5) 

At each point (s, r, (p), there is a unique magnetic surface 

X^x{r,s,(j)) (11.6) 

Inversion of this relation allows us to use now the new coordinate system {s,x,4')- By 
definition of magnetic surface 

B.Vx = (II.7) 

Implies that B has no component in the Vx direction. Chui and Moffatt define a shape- 
function along the magnetic flux tube as R{s, x, 4>) which allows us to define the tube 



geometrically as 

x{s) = X + i?(cos en sin Ob) (II.8) 
Chui-Moffatt used the following fundamental differential form 

dx{s) — eids + e2dx + e^dcf) (II-9) 

Where tie vector basis frame (ei, 62, 63) is given by [3] 

ei = (1 - K cos e)t + Rs{cos e - Rr* sin e)n+ {R, sin 6 + RT*cose) b (11.10) 

Rxicos en + sin 9b) (11.11) 

eg = {R^cos9 - Rsine)n+ (Rcose + R^sine)b (11.12) 
Where r* — r — wiiere L is the length of the tube. The Riemann non-orthogonal 



metric tensor gij where — 1, 2, 3) is given by 
which in matrix format is given by 



(11.13) 



({1-Rk cos ef + + Rs 

Ry^ Rg 

R^T*"^ + Rs^ 



R-^Rg 



R'r*' + RgR^\ 

RxR<i> 



RyRfhR^ + R(h^ R^ + Rft)^ J 



(11.14) 



In the next section one shall see that the poloidal and toroidal magnetic field given by 

Bt = B^ei (11.15) 

Bp = B^e^ (11.16) 
Are constrained by imposing the unstretching magnetic flux tube condition 



(B.V)v= 
onto the magnetic self-induction equation 

dtB={B.V)v+r]V'^B 
Here B"^ vanishes due to the magnetic surface condition above. 



(11.17) 



(11.18) 



Ill Knot energy in unstretched flux tubes 

Just for comparison one writes here the Germano-Ricca [9] Rieniann metric of orthogonal 
coordinates in magnetic flux tube as used in our previous Riemanian flux tubes papers 
[9] is appended 

df = dr' + [r'd^R^ + K''{s)ds^] (III.19) 

where metric factor K[r,s) = (1 — nrcosO) which is contained also in the Chui-Moffatt 
metric, showing that this metric is rather simpler than the Chui-Moffatt one. Despite this 
backdraw, there is a clear advantage to use the more general knotted Riemannian metric, 
and this is the richness of mathematical and physical information, as one shall see in the 
rest of the paper. Of course, Germano-Ricca metric (11.17) is very useful in the solar loops 
topology. The present Ghui-Moffat knotted covariant metric also helps one to investigate 
reconnection of solar tubes by considering thin tubes. Note that by considering the index 
notation (e^) with {A — 1,2) one may express the unstretching condition above as 

(B.V) V = B^dAiv'ei) = B\'dA{ei) = B'^v'pAiej (III.20) 

Note that to simplify matters one considered that the modulus of the velocity flow is 
constant and P^i is the so-called Ricci rotation coefficients, so-commonly use also in 
general relativity. It is clear that these objects here may be written in terms of the Ricci 
coefficients written in terms of the Prenet frame as 

n.dgt — n.K,n — K — Tnss (III. 21) 

n.d,b = -T = T^,t, (III.22) 
t.dsti = = r,,„ (III.23) 

b.dsll =-T = Fhsn (111.24) 

From these last four expressions one infers that the Frenet-Ricci coefficients above have 
the following symmetries 

^jki — —^ikj (III. 25) 

As in Einstein's general relativistic spacetime, where the vanishing of Ricci rotation coeffi- 
cients, imply that the spacetime is Minkowskian or Hat, here the vanishing of Frenet-Ricci 
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rotation coefficients means that the curvature and torsion of the curve vanishes, meaning 
that the hne or magnetic Gelds are Bat. A similar formalism to the one used in this paper 
was used by Ricca [1] to investigate filaments and solitons. From the above expressions 
one may use the following Ricci-like relation 

dA(ei) = r^Aiej (III.26) 

Which clearly shows that 

= (III.27) 

is equivalent to unstretching. One also notes from here, that not all Ricci rotation like 
coefficients F vanish, so one may say that of course, the space does not need necessarily 
be Rat. Next one shall consider, the specific case of Chui-Moffatt covariant metric of 
knotted flux tubes. In the case of the helical knotted tube, the torsion tq and curvature 
kq are constants, a long but straightforward computation leads to the following stretching 
expressions in terms of the Ricci rotation coefficients as 

n.{B.V)v = B\'P Ain.ej (III.28) 

The Ricci rotation coefficients needed can be computed as 
r",i = n.dsei = f,o + (2?^ + t\[Rs + ?^])sine + t\[Rs + ?^]cose (III.29) 

r",2 = n.dse2 = [-R^scose- R^i^^ + r*o] cos^ (111.30) 

^^3 = n.dses = [R^s + - ^*o] cos^ - R4>(^ + ^*o) + Rs) smO (III.31) 

r%2 = t.dse2 = -Ry^K cos 9 (III.32) 

This iast expression has two physically interesting interpretations. Since either, the cur- 
vature K vanishes, which is the trivial straight flux tube, or R^ vanishes, which is the case 
of total constant cross-section tube. Let us now compute the knot magnetic energy of the 
unstretched constant cross-section flux tube, by noting that the RRC could be expressed 
in terms of the Frenet frame directly, instead of the above triad. This allows us to write 

Fnj = ei.dsej = ^iaE^ds[^jplf] (III.33) 



where the vector = {t,n, b) is the compact notation for the Frenet frame. Exphcitly 
this equation yields 

rii2 = ^9,[7ii' + 712' + 713'] (111.34) 

The Ricci coefhcients 'jij can be determined from the values in the expressions for the 
triad frame above. The non-vanishing components of these coefficients are 



7ii = K{r, s) 


(111.35) 


712 = —Rt* sin 9 


(III.36) 


7i3 = Rt* cos 9 


(III.37) 


722 = Rx cos 9 


(III.38) 


723 = Rx ^ 


(III.39) 


732 = —Rsin9 


(III.40) 


733 = Rcos9 


(III.41) 



where in these expressions one already takes into account that R^ and Rg vanish. This 
means that the tube is uniform and its cross-section is constant. Substitution of these 
coefficients into (III. 33) yields 

rii2 = ^[i^:' + T'i?'] (III.42) 

Next one shall see that these affine connections axe very useful to compute the ratio 
between toroidal and poloidal magnetic fields. This can be seen , since according to the 
unstretching condition above 

B^'TiAj = (III.43) 

is equivalent to the equations 

BTaj + B^Tisj = (III.44) 

which in the first case reads 

B^Fiu + B^riS2 = (III.45) 
Now computation of the RRC ri32 as 

ri32 = RRxT* (111.46) 
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yields the following expression for the ratio between toroidal and poloidal magnetic field 
as 

1^ = -6 « -2i?i?,r* (III.47) 

-D J- 132 

Note the this expression shows us that, in the next section example, where the effective 
torsion r* vanishes, in this unstretching case , dynamo action would not be present since 
the toroidal component, would vanish. But note that in the next section the flux tube 
is stretching, so as one shall show the dynamo action can be supported in the knotted flux 
tube. Here one used the approximation of this tubes and helical weak effective torsion 
T* . These expressions allows us immeadiatly to write the expression for the knot energy 
in the Riemannian flux tube manifold[9]. Since the tube is unstretched one may consider 
that its length L is constant. Thus 

M = ^ y v^bii6' + - 2bgi3]dsdxd4> (III.48) 

Since R — R{x) here, the magnetic energy of the knotted MFT becomes 

<M>= -e^L'^Vix) < {B^f > (III.49) 

where e = (^^)^ and X = ffcre Vr is the volume inside the whole tube while is 
the volume under the magnetic surface x — constant. Actually < M > is an averaged or 
mean field energy, and we assume that the average value of B^ squared is constant over 
an averaged volume. Here also 

V{x) = / Vgdsdxd<l> (III.50) 

with the above approximations we were left also with the fact that gis component vanishes 
which also facilitates our computations. Since at the magnetic surface the volume is 
constant from the knot energy, one must conclude that this leads to a marginal dynamo 
action with a constant energy. 
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IV Dynamo action in flux tubes with circular cross- 
section 



The first matematically analyzed fast dynamo model, as a smooth flow of uniform stretch- 
ing, was the suspension of the cat map by Arnold, Zeldovich, Ruzmaikin, and Sokoloff 
[5] in 1981. More recently, Oseledets [13, 14] has been investigated these dynamos in 
random smooth maps, following Baxendale and Rozovskii [15] work. More recently, an 
interesting example of a two-dimensional turbulent dynamo on a Moebius strip flow was 
given by Shukurov, Stepanov and Sokoloff [16], with applications in liquid sodium torus 
dynamo experiment. Chicane, Latushkin and Smith [1 7] have provide a theorem, where 
fast dynamo action could be only be supported in a two-dimensional surface, if this pos- 
sess a constant negative Riemannian curvature. In this paper another new application 
of this theorem, is presented by showing that a fast dynamo action, can be supported 
in inflexionary [18] twisted magnetic flux tube surface of in regions of negative curva- 
ture. This example is similar to chaotic dynamos obtained previously by Boozer [19] 
and Thiffeault and Boozer [20] who investigate eigenvalues, given by kinematic magnetic 
dynamo Lyapunov exponents, instead of Lefschetz numbers. More recently Vishik [21] a 
have presented and proved a anti-fast dynamo theorem, which could be considered as a 
theorem for positive two-dimensional curvature as well. Recently Garcia de Andradc [11] 
has applied Vishik's theorem in dynamo plasmas. This study is important not only to 
kinematic dynamo, but also to hydromagnetic dynamos [22] where Lorentz force, back- 
reaction helps stretching dynamos by plasma flows. In this paper, to simplify matters, 
the resonant condition between poloidal and toroidal frequencies. This condition could 
be tested experimentally in the Perm russian dynamo torus , where the frequencies pro- 
hies could be adjusted. Riemannian negative curvature of Anosov [23] type appears here, 
naturally in inflexionary flux tube surfaces, since the scalar principal toroidal curvature of 
the constant cross-section, is positive, while the curvature in the toroidal eigendirection is 
negative in the deeper regions, where inflexion of the tube takes place. Three-dimensional 
conformal dynamos, [24] have been obtained earlier by stretching of ideal plasma in flux 
tubes. But due to Zeldovich et al [5] and Cowling [10] anti-dynamo theorems for special 
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symmetries of the Qow, two-dimensional dynamo action in incompressible Bows as dealt 
here, are more complex, and deserve further investigation. Here resistivity e does not 
vanish over the two-dimensional Riemannian surface of the flux tubes. The importance 
of investigating the geometry and topology of solar flux tube dynamos, as here, stems 
from the work of M Schuessler [25] on magnetic flux tube dynamos in solar and stellar 
plasma physics. Recently Ashgari-Targhi and Berger [26] have investigated in detail the 
topology of stretch-twist and fold dynamos [27] in flux tubes as well. Let us now consider 
the above CM metric in the case of the circular cross-section, not necessarily constant. 
In this case R(j) vanishes and the effective torsion is also assumed to vanish, which simpli- 
fies computations. Another assumption is that the modulus of the plasma flow velocity 
coincides as and are constants. From the expressions of the last section, one is left 
as 



/(I 


— rk COS 


R^Rg 


\ 











V 








R^j 



(IV.51) 

where now 

ei ^ {1 -K cos e)t + Rser (IV.52) 

63 = R[- sin en + cos 9]b^ Rsq (IV.53) 

where now one has used the relation between the non-orthogonal (ei, 62, 63) and orthog- 
onal basis {t,er,eg) given by 

ee = - sin en + cos e (IV.54) 

and 

Br — cos en + sin e (IV.55) 

Actually, the assumption of circular cross section, could be added the assumption that 
the tube was axially uniform, but since an excess of symmetries can spoil dynamo action, 
here to simplify matters one keep up with the simplest non-orthogonal metric. When 
both assumptions are kept, not only Rg but also Rfj, vanishes. One must note that now 
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the radial component of the magnetic held vanishes. In this case the hux tube basis 
is orthogonal. Besides being stretched by the plasma how, the tubes one has assumed 
above that one has a vanishing effective torsion t* , implies that 

ro='-^ (IV.66) 

which, of course does not mean that the torsion of the tube vanishes, but that it is 
constant. Note also that, now torsion depends directly upon the Gauss linking topological 
number. From the above relation between 9 and s, one is able to obtain the following 
relation between derivative operators 

ro-'ds - -^de (IV.57) 

where ko — ^ is the internal curvature of the tube, which here is positive, since is circular. 
Here Tq is the external radius. The expression kr — ^, represents the internal curvature of 
the tube. An easy computation shows that the relation between the poloidal and toroidal 

frequencies is 

toe = -roBu, (IV.58) 
Note that the resonance hypothesis loq — lOs — <^o implies that 

TO = (IV.59) 

which tremendously simplify the following computations of the chaotic dynamo equation 

dtB = Vx{vxB) (IV.60) 

where diffusive Laplacian term has been suppressed by assuming that, magnetic Reynolds 
number is very high like in ideal astrophysical plasmas environment. By making use of 
the orthogonality between covariant and contravariant basis, given by 

e\e, = S^j (IV.61) 

and the relations 

a^ei = d^es (IV.62) 

dsCi = RssSr (IV.63) 
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and 

dgBs — RgSe — KoRsinOt (IV.64) 
the first term on the LHS of the induction equation reads 

dtB = AB + v^[RB^Rsser + B^{Rsee - RKQsinOt)] (IV.65) 

where one has used the fact that here, the curvature and torsion are constants and coin- 
cident {kq = To), as in hehcal tubes. Since the RHS of the kinematic dynamo equation 
above can be written as 

Vx(vxB) = (B.V)v- (v.V)B (IV.66) 

for divergence-free vector fields, one is able to express the scalar equations which comes 
from the induction equation as 

[Rss{v^ + v^) + XRs]B' = [B\^ - v^B']Rss (IV.67) 

[Rs + ^R]B^ + dsB^ = (IV.68) 
A(l - Rkocos9)B^ - kqB^ = (IV.69) 

To obtain these equations one assumed for simplicity, that the toroidal and poloidal flows 
are equal in modulus. Considering that, there is an cquipartion, between the components 
of the magnetic field and B^ — B^,what is commonly found in plasma physics, one is able 
to reduce these equations to the following system 

2RssV^ + Ai?, = (IV.70) 

A(l - RKocose) - = (IV. 71) 
which yields the following solutions 

f R(a)ds 

B^ = Boe^^'-'^^^ (IV.72) 

and 

Rs = Aoe~^' (IV.73) 
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Note from last equation, that for dynamo action growth rate (A > 0), the integration con- 
stant has to be negative, otherwise the radius would be negative, thus one must concludes 
that the tube shrinks to a constant value. Mathematically this can be expressed as Rg < 0. 
When the tube reaches a constant radius, being suitable to a flux tube approximation, 
the expression (IV. 72 ) for the magnetic field reads 

= Boe^^'-^^ (IV.74) 

A simple physical inspection of this expression shows that the magnetic field grows in 
time, and when the toroidal component of the plasma flow increases, due to previous 
tube shrinking, the magnetic field grows faster than the grow in distance can slow down 
his action. Inclusion of diffusion and back reaction Lorentz forces in the plasma non-ideal 
flow, would certainly change this scenario of a chaotic dynamo. 

V Conclusions 

With the aid of a mathematical formalism popular in Einstein general relativity, called 
Ricci Rotation Coefficients, one imposes a unstretching of magnetic field lines along the 
knotted tube and obtains constraints which allows us to obtain a dynamo action, when one 
moves from one ergodic magnetic surface to another. For uniformly stretched, constant 
cross-section tubes, one obtains a marginal dynamo action, which generalizes the existence 
of non-fast dynamo plasmas for the knotted Riemannian flux tube manifold in comparison 
with the unkotted one given previously in the literature. The formalism applied to the 
dynamo action here can be also used in the investigation of anti-dynamo theorems [10] 
and further investigation in helicity in dynamo plasmas [26] can be further generalized 
to knotted tubes. More complicate knot energy [12] magnetic integrals maybe obtained 
in Riemannian manifolds for more general knotted MFT. Stretching knotted magnetic 
flux, for uniformly stretched, non-constant cross-section tubes, is shown to lead to a 
simple way to obtain chaotic diffusive-free dynamos on an ideal plasma in non- orthogonal. 
Another example of the use of non-orthogonal coordinates in dynamo theory is on general 
relativistic dynamos in rotating stellar objects like torus around black hole plasmas. A 
detailed investigation of these ideas may appear elsewhere [28]. To the interested reader, 
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may found more details about the applications of non-orthogonal cooordinates in plasmas 
may be found in the book of D'haseleer et al [29]. Finally, after I finish this paper I 
became aware of a very recent paper by Maggione et al [30], where they also address 
the CM Riemann metric in magnetic knots, which are defined in terms of the magnetic 
induction equation in ideal plasmas used above. Though the authors do not address the 
issue of dynamos, since their work is more complete and mathematical than ours, would 
be interesting to extend their work, to allow for the presence of dynamo action in magnetic 
knots. 
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